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Abstract 

The superposition of the quantum rotational motion (tunneling) of the encapsulated SC2C2 
complex with the classical rotational motion of the surrounding Cs4 molecule in a powder crystal 
of Sc2C2@C84 fullerite is investigated by theory. Since the quantum rotor is dragged along by 
the C84 molecule, any detection method which couples to the quantum rotor (in casu the C2 
bond of the SC2C2 complex) also probes the thermally excited classical motion (uniaxial rotational 
diffusion and stochastic meroaxial jumps) of the surrounding fullerene. The dynamic rotation- 
rotation response functions in frequency space are obtained as convolutions of quantum and classical 
dynamic correlation functions. The corresponding Raman scattering laws are derived, the overall 
shape of the spectra and the width of the resonance lines are studied as functions of temperature. 
The results of the theory are confronted with experimental low-frequency Raman spectra on powder 
crystals of Sc 2 C2@C 84 [M. Krause et al., Phys. Rev. Lett. 93, 137403 (2004)]. The agreement of 
theory with experiment is very satisfactory in a broad temperature range. 



1 



I. INTRODUCTION 



The existence of an endohedral fullerene, i.e. one or several atoms encapsulated in 
a fullerene molecule, was originally inferred from an analysis of mass spectra of LaCl3- 
impregnated graphite and did lead to the proposal La@C 60 At present the study of 
endohedral metallofullerenes M x @C n , x = 1,2,3,4 and n = 66, 68, 72, 74, ... , 100, where M 
are group II and III metals such as Sc, Y, ... or lanthanides Ce, . . ., Lu, is a subject of 
interdisciplinary research [2] in physics, chemistry and materials sciences. By now one is 
able to produce materials where not only single atoms but clusters of atoms are encapsu- 
lated [3(. Due to charge transfer between the cluster and the surrounding carbon cage it is 
possible to obtain molecular-like complexes which do not exist otherwise (i.e. in absence of 
encapsulation) and which have unusual properties. Not only clusters with metal atoms of a 
same kind, such as the Sc 3 trimer in Sc 3 @C 82 are produced, but also clusters composed of 
different kinds of atoms. A remarkable case is the production of Sc2C2@C§4 in crystalline 
powder form |4j. The powder crystal is composed of crystallites where the Sc2C2@Cs4 units 
are arranged with average space group symmetry Fm3m. From spectroscopic and struc- 
tural characterization by NMR- and synchrotron X-ray diffraction experiments [J] it follows 
that the Sc 2 C 2 complex is encaged as a rigid unit with point group symmetry D 2 h in a C 8 4 
fullerene of symmetry D 2 d (isomer III, number 23 [5|). The center of mass of SC2C2 coin- 
cides with the center of mass of the molecule. The two Sc atoms are located at a distance 
of 4.29(2) A on the long C2 (S4) axis of C 8 4. The two C atoms of SC2C2 are located in the 
plane containing the two C2 axes perpendicular to the long axis of the C 8 4 molecule and 
have a calculated distance of 1.28 A. This distance lies between those of typical double and 
triple carbon bonds, and is consistent with the experimental and calculated C-C stretching 
frequency of the C2 unit (exp. 1745 cm -1 , calc. 1742 cm -1 ) [6!]. In the following we will 
speak of this C-C bond as a C2 unit or molecule. Indeed low energy Raman spectra 6| on 
powder samples of Sc2C2@C 8 4 in a temperature range 25 - 150 K (Kelvin) have revealed the 
existence of quantized rotational states of the C2 unit. The Raman lines' positions reflect 
transitions between energy levels of a C2 planar quantum rotor in a fourfold static potential 
due to the surrounding C 84 cage. Therefore, one can speak of a quantum gyroscope. 

In Ref. [^I the potential parameters of the encaged quantum rotor were obtained from 
density functional calculations using the VASP (Vienna ab initio simulation package) code 
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The energy levels were then determined from the solution of the Schrodinger equation. 
Within this approach, the Raman spectra consist of infinitely sharp lines while experimen- 
tally the lines are broadened and have a characteristic temperature behavior. A reason for 
this shortcoming is the restriction of the role of the encapsulating C§4 molecule to a purely 
static body. Since the measured transition frequencies are in the range of 10 - 80 cm -1 and 
since the line broadening is of the order of a few cm -1 , any involvement of internal vibra- 
tional modes of the Cs4 cage as well as of stretching modes Sc-Cs4 can be excluded. The 
latter are of higher frequencies and have been measured in C 8 4 and in Sc 2 @C 8 4 by infrared 
and Raman techniques . However the low-frequency external rotational modes of the 
C§4 molecule and their superposition with the transitions of the quantum rotor should be 
retained: indeed the encapsulated Sc 2 C 2 gyroscope is dragged by the classical rotational 
motion of the Cs4 molecule and this dragging will affect the Raman spectrum of the C 2 
unit. It follows that an experimental probe such as Raman scattering which couples to the 
encapsulated species in an endohedral complex, in casu Sc 2 C 2 , also yields information on 
the dynamics of the encapsulating molecule, in casu Cs4- 

In the present paper we will extend the theoretical interpretation given in Ref. O and 
develop a unified theory where the quantum mechanical motion of the Sc 2 C 2 complex is cou- 
pled to the thermally excited classical rotational motion of the Cs4 fullerene. The coupling 
results from the fact that the long axis of the quantum gyroscope coincides with the £4 axis of 
the surrounding Cs4 molecule. In a given crystallite the C 8 4 molecules are randomly oriented 
with their long C 2 axis in equivalent (100) directions of the face-centered cubic (fee) unit cell 
4j. We call this situation meroaxial disorder (this terminology seems to be more appropriate 
than merohedral disorder). We start from a model where at low temperature the meroax- 
ially oriented C 8 4 molecules in the fee crystal perform uniaxial rotational diffusions about 
their long axis. Such a classical motion can be seen as a time-dependent modulation of the 
fourfold potential experienced by the quantum rotor and causes a temperature-dependent 
broadening of the quantum levels. An additional broadening effect is to be expected from 
the stochastic reorientations of the C 8 4 molecules among the meroaxial directions which 
should become increasingly important at higher temperature. In addition, the stochastic 
reorientations lead to the appearance of a temperature-dependent quasi-elastic peak in the 
Raman spectrum. 

The content of the paper is as follows. In Section II we write down the Raman scatter- 
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ing law for the C2-unit belonging to the Sc 2 C2 complex encapsulated by the C 8 4 molecule 
in Sc 2 C 2 @C 84 fullerite. We start from a single crystal with Fm3m structure and static 
meroaxial disorder of the C 8 4 molecules. Assuming a quantum mechanical rotational mo- 
tion of the SC2C2 complexes and a classical rotational motion of the C 8 4 molecules, the 
dynamic polarizability-polarizability correlation function is decoupled in a product of corre- 
lation functions for the rotational dynamics of Sc 2 C 2 and C 84 respectively. The scattering 
law is obtained as a convolution of these correlation functions in Fourier space. Next (Sect. 
Ill) we calculate the correlation functions, using a quantum mechanical tunneling model 
for the Sc 2 C 2 complex and a uniaxial rotational diffusion model for the surrounding C 8 4 
molecule. The rotational diffusion motion of the encapsulating C 8 4 molecule leads to a lin- 
ear temperature-dependent broadening of the energy transition lines of the C 2 planar rotor. 
In Sect. IV we extend the theory to a powder crystal which consists of arbitrarily oriented 
crystallites, each with Fm3m space group symmetry and static meroaxial disorder. In the 
following (Sect. V) we consider the case of dynamic meroaxial disorder, describing the re- 
orientations of C 8 4 molecules among the three meroaxial directions by a stochastic jump 
model. This model yields an exponential temperature-dependent broadening of the transi- 
tion lines. In the last Section VI we give a numerical evaluation of the Raman scattering law 
for a Sc 2 C 2 @C 8 4 powder crystal where quantum mechanical tunneling of the encapsulated 
Sc 2 C 2 units is superimposed by uniaxial rotational diffusion and dynamic meroaxial disorder 
of the C 8 4 molecules. The temperature dependence of the line intensities and of the line 
broadenings is discussed. 



II. RAMAN SCATTERING LAW 



We will derive the Raman scattering law where we limit ourselves to the interaction of 
the incident laser light with the plane rotational motion of the induced dipole of the C-C 
bond belonging to the Sc 2 C 2 complex of Sc 2 C 2 @C 84 . This means that we consider the low 
frequency part of the spectrum (say < 100 cur 1 ). The Sc 2 C 2 complex is centered in the 
origin (center-of-mass position) of C 84 . The long axis of Sc 2 C 2 coincides with the £4 axis of 
C 8 4. The C-C bond of Sc 2 C 2 lies in the plane containing the secondary C 2 axes of C 8 4 and 
rotates about the S4 axis. In that respect we will consider the C-C bond as a C 2 planar 
rotor which experiences a fourfold potential inside the C 8 4 molecule. Our formulation of 
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FIG. 1: View of the Cs4 molecule along the S4 axis. 



the Raman scattering law is an extension of the conventional theory [10|, in as much as 
we describe a situation where the rotational motion of the induced dipole with respect to 
the laboratory-fixed frame is a superposition of the quantum motion of the C2 planar rotor 
inside the Cg4 molecule and of the classical motion of the C$4 molecule in the laboratory 
frame. 

We start with considering a single crystal of Sc 2 C 2 @C84 units with static meroaxial dis- 
order. We assume that the Sc 2 C 2 @C 84 units are statistically independent, hence it will be 
sufficient to consider one single representative unit. The cubic crystal axes (X', Y', Z') are 
chosen to coincide with the laboratory- fixed cubic coordinate system (X, Y, Z) . We consider 
a cubic system of axes (£, 77, () fixed in the Cg4 molecule such that the £ axis coincides with 
the 5*4 axis while 7? and ( coincide with the secondary twofold axes (Fig. [T]). The meroaxial 
orientations of the Cs4 molecules correspond to the situation where the £ axes are randomly 
oriented along the X', Y' or Z' crystal axes. The C 2 units then rotate in the planes (Y 1 , Z'), 
(Z',X f ) or (X',Y f ) respectively (Fig. [2]). We say that the C 8 4 molecule is in standard ori- 
entation if the S4 axis coincides with the laboratory-fixed X axis and the plane containing 
the secondary C 2 axes coincides with the laboratory (Y, Z) plane. The ( axis forms an angle 
v with the Z axis, while the C-C bond forms an angle r with the ( axis. Hence the polar 
angle 9 of the C-C bond with the laboratory Z axis (Fig. EJ) is a sum of two terms: 

e = u + T. (2.1) 

Since the C 2 rotor is confined to the (Y, Z) plane, the azimuthal angle <fi measured away 
from X has value n/2. The distinction of two contributions to the angle 9 is essential. 
In the following we will assign the angular variable r to the quantum mechanical tunneling 
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FIG. 2: The SC2C2 complex in the crystal-fixed cubic coordinate system (X' , Y', Z') while the 
molecule is in standard orientation. 
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FIG. 3: Orientation of the C2 bond of SC2C2 in the rotatory reflection plane of the Cg4 molecule 
(C84 in standard orientation). 



of the SC2C2 complex about its long axis inside the Cs4 cage and the angular variable v to the 
thermally excited classical rotation of the C§4 molecule about the 64 axis. The assumption of 
classical uniaxial rotational diffusion motion as a first approximation to the dynamics of the 
C»4 molecule at low temperature is motivated by the structural results of meroaxial disorder 
|4|]. It is also inspired from the dynamics of solid C70 in the rhombohedral and monoclinic 
phases. There the importance of uniaxial rotational diffusion about the long axis of the C70 
molecule has been probed by muon spin spectroscopy 
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131 ]. nuclear magnetic resonance 
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UM, lly] and inelastic neutron scattering [17]. 



We treat the C-C bond of SC2C2 as a rigid cylindrical rod with longitudinal and transverse 
static polarizability an and a± respectively. The Raman scattering law for incident and 
scatterd radiation in Z direction is given by the Fourier transform of the time- dependent 
autocorrelation function of the polarizability a zz : 

Rzzzz(uj) = — / dte^(a zz (t)a zz (0)). (2.2) 

Here N is the number of Sc2C2@Cg4 units and u is the frequency difference of incident 
and scattered radiation. The polarizability has to be understood as an average over the 
three meroaxial molecular orientations. In the following we label these orientations by a 
superscript (i), i = 1,2,3. If the Cs4 molecule is in standard orientation (£ axis || X), or 
in orientation £ || Y, the corresponding orientation-dependent polarizabilities a zz and a zz 

4l = a f z = ^l + ^ ws2e , (2 .3) 

while with £ || Z one has 

af z = a ± , (2.4) 

independent of 9. In the case of meroaxial disorder, the average polarizability is given by 

3 



1 v — -\ 2 

"zz = -^«zz = a + 3 &cos261 ( 2 - 5 ) 
i=i 

where we have defined 



Oi\\+2a ± Q!||-Q!_L 

a = ^^' b = ^T~ ■ (2 ' 6) 

Hence the time-dependent correlation function reads 

Ab 2 

(a zz {t)a zz {0)) = a 2 + — (cos 29 (t) cos 29(0)). (2.7) 

y 

Similarly we obtain for incident radiation in Z direction and scattered radiation in Y 
direction 

Rzyzy^) = — / dte^(a ZY (t)a ZY (0)). (2.8) 
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If the C84 molecule is in standard orientation, 

a ( z ) Y = bsm29, (2.9) 

while for £ || Y, af], = 0, and £ || Z, a^ Y — 0- The average polarizability for the case of 
meroaxial disorder reads 

a Z Y = ^sm26, (2.10) 

and the correlation function becomes 

b 2 

(a ZY (t)a ZY (0)) = -(sin 20(i) sin 20(0)). (2.11) 

y 

The problem of determining the scattering laws Rzzzzij^) an d Rzyzy(u) consists in the 
calculation of the orientation-orientation thermal correlation functions 

C(t) = (cos 26 (t) cos 26(0)), (2.12) 

S(t) = (sin20(*) sin 26(0)}, (2.13) 

Taking into account the basic relation Eq. (12. ip . we expand in terms of cos2r, sin2r, cos2z/ 
and sin 2v thereby obtaining correlation functions of the form 

C cccc (t) = (cos2r(£) cos2z/(t) cos2r(0) cos2z/(0)), (2.14) 

S cscs (t) = (cos2r(t)sin2z/(t)cos2r(0)sin2z/(0)), (2.15) 

and similarly for C ssss (t) and S scsc (t). Observing that r refers to quantum dynamics of C2 
and v to classical dynamics of C84, we decouple the thermal averages: 

C cccc = Q cc (t)F cc (t), C ssss = Q ss (t)F ss (t), (2.16) 

S cscs = Q cc (t)F ss (t), S scsc = Q ss (t)F cc (t). (2.17) 

Here the correlation functions 

Q cc (t) = (cos2r(t) cos2r(0)), (2.18) 

Q ss (t) = (sin2r(t)sin2r(0)), (2.19) 

describe the quantum dynamics of the C2 unit while the correlation functions 

F cc (t) = (cos2u(t) cos2z/(0)), (2.20) 

F ss (t) = (sin2z/(t)sin2z/(0)), (2.21) 
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describe the classical dynamics of the surrounding Cs4 molecule. Finally quantum and 
classical dynamics occur as products of correlation functions: 

C(t) = Q cc (t)F cc (t) + Q ss (t)F ss (t), (2.22) 

S(t) = Q ss {t)F cc {t) + Q cc {t)F ss {t). (2.23) 

Defining Fourier transforms 

i n+oo 

Q(u) = — / dte^Qit), (2.24) 



2tt 



— oo 
oo 



F(u) = ^-I dte^F(t), (2.25) 

and using Eqs. (I2.22p . (I2.23p . (12. 7p and (12.111) . we rewrite the Raman scattering law in terms 
of convolutions of Fourier-transformed quantum and classical correlation functions, thereby 
obtaining 

Ab 2 

Rzzzz{uj) = N[a 2 8(u) + -^-C(u)], (2.26) 



with scattering function 



and 



/+oo 
duo 1 [Q cc {u - uj')F cc (uj') + Q ss {u - u/)F"(u/)] , (2.27) 
-oo 



b 2 

Rzyzy{uj) = N-S(u), (2.28) 
y 



with scattering function 

/+oo 
dJ [Q ss {u - uj')F cc {uj') + Q cc {uj - u')F ss {u')} . (2.29) 
-oo 

The first term in brackets on the right-hand side of Eq. (12.261) corresponds to the unshifted 
Rayleigh line of the spectrum while the function C{uj) (as also S(uS) in Eq. ( 12.28(1 ) accounts 
for the inelastic part. Expressions (12. 27(1 and (12.291) which are convolutions in Fourier space 
show that the quantum motion of the C2 rotor is modulated by the classical rotational 
motion of the surrounding C§4 cage. This is an example of "direct coupling" of two motions 
through the detection process [18j] , in contradistinction to the "indirect coupling" through 
a Hamiltonian. The origin of the direct coupling here is due to the fact that the detection 
angle 9 is a sum of two terms, Eq. (12. ip . 

In the next section we will calculate the quantum mechanical and classical orientational 
correlation functions for C 2 and C 84 respectively. 



III. DYNAMIC CORRELATIONS 



A. C2 quantum rotor 

The quantum mechanics of a diatomic molecular rotor in crystals goes back to Pauling 
[3]. A still valid review of the subject of single particle rotations in molecular crystals has 
been given by W. Press 20] . We will calculate the orientational autocorrelation functions 



burfold 



20j,|21|. 



Q cc and Q ss by starting from the model of the C2 planar quantum rotor in the 
potential due to the C§4 cage. We will refer to this motion as rotational tunneling 
We will show that the resonances of the correlation functions Q cc (u>) and Q ss (u) are due 
to transitions between tunneling energy levels. The sole degree of freedom is the angle r 
which accounts for the rotatory motion of C2 with respect to the cage. The corresponding 
Schrodinger equation reads 



-£-^t + ^(1-4cos4t) 

CLT Z 2 



%1)(t) = Eip(r). 



(3.1) 



Here B = h 2 /2I is the rotational constant and / the moment of inertia of C2, Vq is the 
barrier height of the potential. The rotational constant has the dimension of an energy, 
from experiment [6j one deduces B = 1.73 cm^ 1 (wave number units) and Vo = 8B. These 
values are supported by ab initio density functional calculations 6]. Equation ( 13.1 ft which 
is an extension of Mathieu's equation 19|, [22] is also called Hill's equation [23]. With the 
definitions 

V \ Vo 



1 

B 



E- 



AB'' 



Eq. (EUD reads 



dr 2 



+ a + 2gcos4r 



V>(t) = 0. 



(3.2) 



(3.3) 



From symmetry considerations (nuclear spin is zero for 12 C, electron wave function of C 2 , - 
is totally symmetric), it follows that the rotational wave function ip(r) must be symmetric 
with respect to the operation r — > t — n. For even periodic solutions one makes the ansatz 



^ + ( r ) = ^ A 2m cos(2mr), 



(3.4) 



m=0 

m — 0, 1, 2, Equation (13.31) then leads to an infinite system of homogeneous equations 

for the coefficients v4 2m . Truncation of this system for a given value m = N leads to N + 1 
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equations which separate into two systems: a first one for { Aq, A±, . . . , A2n} arid a second 
one for [A 2 , A e , . . . , A2N-2} (we take N even). Solving for the two discriminants yields the 
roots ot2 m (q) for m — 0, 2, . . . , N and m — 1,3, ... , N—l. In case of zero potential, i.e. q — 0, 
these solutions reduce to the free planar rotor energies {-Egm^ = 0)} = {0, . . . , (2m) 2 B, . . .} 
with normalized eigenfunctions 

in the interval < r < 7T. The ansatz for odd periodic solutions reads 

00 

if)~{r) = 2J -B2mSin(2mr). (3.6) 

m=l 

Proceeding as before one determines the roots a^" m (g). In case of zero potential the eigen- 
functions are 

r , / \ ~i |sin2r sin(2mr) I . 
{ ^. W } = J-^-,...,_^2,...j. (3.7) 

In the following we will label the energy eigenfunctions and eigenvalues by the double index 
(cr, 2m), er = ±, also in the case of nonzero potential. In Fig. H] we show plots of E2 ™^ = 
a.l m {q) + 2q. 

We next perform a spectral decomposition of the correlation functions Q cc {t) and Q ss {t) 
in terms of eigenfunctions and eigenvalues of the Schrodinger equation (13. ip . In general form 
the result reads 



where 

Cii = (i|cos2rU), (3.10) 
5 y = (i|siii2r|j). (3.11) 

Here the label z (j) stands for the double index (a, 2m) of the solutions of the Schrodinger 
equation. We calculate the matrix elements CV,- and Sij with the free planar rotor energies. 
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FIG. 4: Energy levels of the C2 planar quantum rotor in the fourfold molecular potential as a 
function of the potential strength q (dimensionless units). Only levels up to 2m = 8 are shown. 



Symmetry implies that only functions of a same parity (+, +) or 



contribute to C 



while Sij differs from zero only for functions i, j with different parity. For instance 



kJ 2m2n 



w cos(2mr) cos(2nr) 1 
dr , — _ cos2r ^==- = -d m>n ±i, 



V^/2 2 



cos(2mr) sin(2nr) 1 
dr : — _ sm 2r 1 ==- = -5 m>n±1 , 



(3.12) 
(3.13) 



These matrix elements imply selection rules for transitions between energy levels. We take 
Fourier transforms of Eqs. (13. 8p and (13.91) . using the identity 



— dte^e^- E ^ h = d(uj-( ^— ^ 



(3.14) 



We insert the energies 



E° m = Ba° m {q) + 



(3.15) 



and take into account the selection rules (13.121) - (I3.13p . Defining the frequency transfer 

(3.16) 



T?0 T?a 

- Cj 2n Cj 2m 



h 



12 



we obtain 



i f 00 — E + IT 

m=l 



+ 7) [5(^-^12 ) + Z^ 9 - u mm+l) + S ( u - u mm-l)\ f> ( 3 - 17 ) 



2 

m=2 



with 



00 

Z = e~ E o/ T + (e" E -/ T + e- £ -/ T ) . (3.18) 



?Ti=l 



Similarly we get 



2 

m=2 



^ - P 2™' r -1 



m=l 

We notice that in absence of the uniaxial rotation of the Cs4 cage, i.e. for v = 0, the 
correlation functions Eqs. (I2.20p and ( 12. 21ft reduce to constants: F cc = 1, F ss = 0. Hence 
the spectral functions C(uj) and S(uj) entering the Raman scattering laws Eqs. (12.261) and 
( I2.28P reduce to Q cc (uj) and Q ss (u>) and exhibit infinitely sharp 5-peaks which account for 
transitions between quantized planar rotor states. In Table I we have quoted some values 
(m < 4) of w™ n for a = (free rotor) and q = 2 (value of the potential strength taken from 
experiment in Ref . |6|) . 



B. uniaxial rotational diffusion 

In order to calculate the classical correlation functions F cc (t) and F ss (t) we treat the C§4 
molecule as a classical uniaxially diffusing rotor with rotation axis 5*4 in coincidence with 
a cubic crystal axis, in casu the X' axis. The corresponding rotation angle v is measured 
away from the Z' axis. Equivalently one considers the S4 axis along Y' and Z' (meroaxial 
disorder). Given the S4 axis it would be tempting to study this motion in a crystal field 



potential of fourfold symmetry. Such a study can be carried out along the lines of Ref. 



24 



and 



leads to a continued fraction expansion in terms of frequency moments of the orientational 
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TABLE I: Tunneling frequency transfers u}™ n , n = m db 1, a, a' = ±, in units cm 



111 


UJ ++ 
vnvn+1 








mm—1 


^mm— 1 


mm—l 




q = 





6.92 




6.92 












1 


20.76 


20.76 


20.76 


20.76 


-6.92 






-6.92 


2 


34.60 


34.60 


34.60 


34.60 


-20.76 


-20.76 


-20.76 


-20.76 


3 


48.44 


48.44 


48.44 


48.44 


-34.60 


-34.60 


-34.60 


-34.60 


4 


62.28 


62.28 


62.28 


62.28 


-48.44 


-48.44 


-48.44 


-48.44 


q = 2 





4.10 




10.99 












1 


25.12 


17.39 


24.28 


18.23 


-4.10 






-10.99 


2 


34.00 


34.87 


34.02 


34.84 


-25.12 


-17.39 


-18.23 


-24.28 


3 


48.40 


48.38 


48.40 


48.38 


-34.00 


-34.87 


-34.84 


-34.02 


4 


62.26 


62.26 


62.26 


62.26 


-48.40 


-48.38 


-48.38 


-48.40 



variables. It is adequate in the case of a strong crystal field potential since then one can limit 
the continued fraction to a few steps. However this approximation is not valid in the case 
of weak potentials. Since the equator of the Cs4 molecule for rotations about S4 deviates 
only slightly from circular shape, we prefer to consider the rotator about the S4 axis in 
the rotational- diffusion approximation. This model has the obvious advantage of simplicity 
and leads to a linear temperature-dependent broadening of the tunneling transition lines. 
Within the uniaxial diffusion model the Cs4 molecule experiences a random rotational torque 
(also called Brownian motion torque) about its S4 axis. This torque is caused by the thermal 
motion of the surrounding lattice (heat bath). In that respect the present problem is different 
from the situation of the heavy symmetrical top with gravitational torque since on the 
molecular scale the effect of gravitation is negligible in comparison with the heat bath. 

The idea of rotational diffusion goes back to Debye [25[] who applied the co ncep t of 
rotational Brownian motion to the theory of dielectric relaxation (see also 2o| and 26]). In 
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Appendix [A] we give some details for the present problem. As a result we obtain 

F cc (t) = (cos2i/(t) cos2z/(0)) = ]^r^\ 
F ss {t) = (sin2i/(*)sin2i/(0)) = i e ~ 4I?R *. 

Here the rotational diffusion coefficient Dr is given by the Einstein relation 



D 



R 



c 



(3.20) 
(3.21) 

(3.22) 



where ( is the friction coefficient and T the temperature. The equality of F cc (t) and F ss (t) 
is a consequence of our neglect of the crystal field potential within the large-friction approx- 
imation. From Eqs. ( 12. 22ft and ( 12.231) one sees that then 

C(t) = S(t). (3.23) 

In the following we will neglect the superscripts ss and cc on F cc and F ss and write just F . 
The Fourier transform is obtained from Eq. (13.211) with the result 



Zn 



AD 



R 



u 2 + 16Dl_ 

We rewrite the scattering functions Eqs. (I2.27P and (12.291) as 

/+oo 
dJ [Q ss {uj - J) + Q cc {u - uj')]F{uj'). 
-oo 

Using expressions (13.171) and (I3.19P we carry out the integral over to' and obtain 

C{u) = C ++ (uj) + C__(o;) + C + 4uj) + C_ + (uj) 

where 



(3.24) 



(3.25) 



(3.26) 
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-) + E 



[F{u - ^+++i) + F(u - utL-i 



m=l 



(3.27) 
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1 e"^ 



2Z 



1 



V1 e -^ /T 



[F(u - uj mm+1 ) + F(^j - u mm _ x 



m=2 



e - S 2 + / T 



(3.28) 
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E 



2Z 1 ^ 2 

. m=l 
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mm— 1 



(3.29) 
(3.30) 
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We see that C(u) is a sum of weighted Lorentzians 

F(u - u™ m±1 ) = — 
lix 

centered around the allowed frequency transfers u = tu^^i an d of width 8Dr (full width 
half maximum). Since Dr has dimension s" 1 , it follows from Eq. ( 13. 22ft that £ has the 
dimension of an action. We write ( = ( n h, where ( n is a dimensionless number taken as 
parameter. We then obtain Dr = 2.08 x 10 10 (T/( n ) s _1 where T is measured in Kelvin. 
Equivalently, Dr = 0.694(T/£ n ) cm -1 . Since to our knowledge there are so far no direct 
measurements of the orientational dynamics of the C§4 molecule in Sc2C2@Cg4, we will choose 
a value of Dr such that the correlation time r c = (4Dr) _1 has a value that is intermediate 
between the values of 2 ns and 5 ps measured by NMR experiments for the C70 molecule in 
;he low-temperature monoclinic and high-temperature fee phases of solid C70, respectively 
1^ | . Assuming that ( n = 100 is a realistic value (then Dr = 10 10 s _1 at T = 50 K), we 
have plotted the scattering function C(u) for several temperatures in Fig. [5j The resonances 
are centered at the frequency transfers uj™^ for the potential strength q = 2. The spectra 
reflect the characteristic assymetries for uj > and u> < due to anti-Stokes and Stokes 
processes, respectively. In our calculations, we have included transitions with the values 
m, n — 0, . . . , 19. 

We notice that if one artificially excludes the tunneling motion of the C2 unit by taking 
a fixed value, say 0, for the angle r, one finds that Q ss (uj) = and Q cc = S(u). Then Eq. 
(I3.25P becomes C(u) = F(u). Since the C 2 -unit is dragged along with the classical rotational 
diffusion of the encapsulating Cs4 molecule, its polarizability is changing accordingly with 
time and the Raman scattering laws Rzzzz(u) and Rzyzy{^) will exhibit a Lorentzian 
F(u) of width ADr centered at uj = 0. 



AD 



R 



[UJ 



(3.31) 



IV. POWDER AVERAGES 

In Sect. IHII we have considered a cubic crystal with crystal axes (X', Y', Z') in coincidence 
with the laboratory- fixed cubic axes (X, Y, Z). Since experiments are performed on powder 
samples, we will extend the previous results. The powder sample consists of a large number 
of arbitrarily oriented cubic crystallites, each crystallite has symmetry Fm3m where the 



Sc2C2@Cg4 units are meroaxially disordered 



4J. We first will consider one single crystallite 
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FIG. 5: Scattering function C(uj) for T = 20 K, 60 K and 120 K. The width 8Dr is the same for 
all resonance lines and increases linearly from 1.11 cm -1 to 6.66 cm -1 for T = 20 K and 120 K, 
respectively. The inset at T = 20 K shows the splitting of the and uof^ resonances. 
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where the crystal-fixed system of axes is related to the laboratory- fixed system of axes by 
the Euler angles (a, (3, 7). The C2 rotors are now moving in planes normal to the X', Y', 
Z' axes of the rotated coordinate system. This means that the polarizabilities z or a ZY 
measured in the laboratory-fixed coordinate system will depend on the Euler angles of the 
given crystallite. 

In Appendix |B] we have calculated the polarizability components a zz which are obtained 
from a zz by application of the rotation operation i?(a,/3, 7 ): 

af z R = R(a,P, 1 )af z . (4.1) 



The meroaxial average 



is obtained as 



«fs = a + \Yi[A%(P,i) + Bf z {^l) sin 29 + C® (/3, 7 ) cos 29]. (4.3) 



3 * , 



where % = 1 refers to £ || X f , % = 2 to £ || Y 7 and z = 3 to £ || Z'. The coefficients A^(/3, 7), 
-E>^(/3, 7 ) and C^(/3, 7 ) are derived in Appendix [Bj they are found to depend on only two 
Euler angles. 

In the present section we assume that the meroaxial disorder is static or equivalently 
there are no reorientations of the C§4 molecules among the three meroaxial directions in 
a given crystallite. The angle 9 is then the sole dynamic quantity. The time- dependent 
polarizability correlation function per molecule in the given crystallite is obtained as 

<a£ z (t)a£ z (0)> = « 2 + f E l A zz^ z + B®B<p z S(t) + Cf z Cf z C{t)] . (4.4) 

The correlation functions C(t) and S(t), defined by Eqs. f)2.12p and (12.131) respectively, have 
been calculated in Sects. II and III. 

The powder average for a function F((3,j) is defined as 

1 r^n 







F= — I dy I d(3 sin (3 F ((3, y). (4.5) 



The results for the products A ZZ A ZZ , B ZZ B ZZ , . . . are quoted in Appendix [Bl The powder- 



averaged polarizability correlation function per molecule reads 

(4.6) 



(«fzW«fz(0)> = « 2 + ¥ 
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Taking into account S(t) = C(t), Eq. (13.231) . we obtain the Raman scattering law for a 
powder-averaged sample with meroaxial disorder: 

Rzzzz{uj) = N (a 2 5(co) + ^b 2 C(u)^j . (4.7) 

The expression for a single crystal with meroaxial disorder has been given by Eq. (I2.26p . 
In an analogous way one calculates 

3 



,,R 

3 



with the result 

3 



a 



*y = \ E [ A %r( a > A 7) + Bf Y (a, (3, 7 ) sin2# + C®(a, A 7) cos 20] . (4.9) 



3 * i 



The coefficients A ZY , . . . , C ZY are given in Appendix [Bj The time-dependent polarizability 
correlation function per molecule reads 

{a R ZY {t)a R ZY (0)> = | [ A zy A zy + Bf Y Bf Y S{t) + C« <7&C7(t)] . (4.10) 
The powder average of a function F(a,/3,j) is defined by 

i /*2"7r /*27T /*7r 

F = — da d-f dp sin/?F(a, /?, 7 ). (4.11) 
8tH i Jo 7 



Taking into account the powder averages A ZY A ZY etc., calculated in Appendix [Bj we 
obtain 



(4.12) 



The Raman scattering law then reads 



Rzyzy(uj) = N-C(u), (4.13) 

where again we have used S(t) = C(t), Eq. (13.231) . We see that the powder-averaged 
expression is the same as the one for a single crystal with meroaxial disorder, Eq. (12.281) . 
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V. DYNAMIC MERO AXIAL DISORDER 



So far we have assumed that the orientation of the long axis (S4) of the Cs4 molecule in 
a given cubic crystallite along the equivalent (100) directions is random but static. The sole 
effect of the heat bath was the uniaxial rotational diffusion studied in Sect. IIII Bl While 
this situation of static meroaxial disorder is realistic at temperatures inferior to say 100 
K, it becomes less valid at higher T. Here again we refer to the situation in solid C70 
where with increasing temperature it is found that the uniaxial rotation axis flips between 
different symmetry equivalent orientations such that the rotational motion becomes more 



n 



and more isotropic [l2|, ll3|, llJ, ll5|, Il6, \y\ . We therefore will extend the previous model 
and take into account the situation where a molecule at a given lattice site in one crystallite 
changes orientation with the S4 axis jumping randomly between equivalent potential minima 
in (100) directions. Here the heat bath causes stochastic torques about axes perpendicular 
to the long axis of the C§4 molecule or equivalently perpendicular to the rotation axis of the 
encapsulated SC2C2 quantum gyroscope. We recall that accordingly the normal to the plane 
of the C2 quantum rotor will change its orientation. Within a simple three sites stochastic 



jump model (see e.g. 



20||). the conditional probability p(i,t\j,Q) to find a C 8 4 molecule in 



an orientation % = 1, 2, 3 at time t > when it was in orientation j = 1, 2, 3 at time is 
obtained by solving a system of three linear differential equations. One obtains 

p(i,t\j, 0)= 1 -(l + 2e- 3wt ),t=j, (5.1) 
p(i,t\j,0) = ^(l-e- 3wt ),i^j, (5.2) 

where w is the transition rate for a molecular reorientation. We associate the transition rate 
with the inverse of a relaxation time: 

w = - = -e~ E ^ T . (5.3) 

Here we have assumed an Arrhenius-type law, known from reaction rate theory 2a, |29| . 
where l/r is an attempt frequency and E a an activation energy for meroaxial reorientations 
of the Sc 2 C 2 @C 8 4 complex as a whole. The equilbrium value of the conditional probability 
is independent of the initial and final orientation and corresponds to an a priori probability: 

lim p(i,t\j,0) = l-. (5.4) 

t — >oo J 
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In the previous section the meroaxial orientations within a given crystallite have been char- 
acterized by the coefficients {A§ z , Bf z ,Cf z }, {A%, B { ZY ,C^ Y } in Eqs. (E3D and ([U]) of 
the polarizabilities ot zz and a ZY . Treating these coefficients as dynamic stochastic quantities 
we obtain instead of Eqs. (I4.4p and (I4.10p 

<«fz(*)«Sz(0)> = a 2 + b 2 [(A zz (t)A zz (0)) 

+ (B zz (t)B zz (0))S(t) + (C zz (t)C zz (0))C(t)] , (5.5) 
(a* Y (t)aU0)) = b 2 [(A ZY (t)A ZY (0)) 

+ (B Z y{t)B Z y(0))S{t) + (C Z y(t)C Z y(0))C{t)] . (5.6) 

The correlation functions \A zz (t)A zz (0)), . . . , (C Z y(t)C Z y(0)) which refer to meroaxial 
reorientations are calculated within the frame of the stochastic jump model. For instance 
for a given set [AW, i = 1, 2, 3} one has 

(A(t)A(O)) = ± £ ^A^p{i, t\j, 0), (5.7) 

where the conditional probabilities p(i,t\j,0) are given by Eqs. ( 15. II) and ( 15. 21) . while the 
factor 1/3 accounts for the equilibrium initial probability. Since the coefficients A^ depend 
on the Euler angles which specify the orientation of a given crystallite (Sect. lTVT) . the powder- 
averaged correlation functions are obtained by averaging over the Euler angles: 

(A(t)A(0)} = ^A®A®p{i, t\j, 0). (5.8) 

Taking into account the values of the powder-averaged products given in Appendix [Bj we 
obtain: 



A zz (t)A zz (0)) = -e" 3 ^, (5.9) 



[B zz (t)B zz (0)) = JL [1 + 2e - 3 ^] , (5.10) 
12 



[C zz (t)C zz (0)) = — [1 + 3e- 3 ^] . (5.11) 



The powder average of Eq. ( 15. 5p then reads 



Ab 2 



<af z (t)af z (0)> = a 2 + — D{t), (5.12) 
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where the function D(t) is given by 

Q 1 Q 

D(t) = [C(t) + -e^ t/T + — C(t)e- 3t/T ] . (5.13) 



Here we have used again C(t) = S(t), Eq. (13.231) . The first term on the right-hand side C(t) 
accounts for the superposition of the quantum motion (tunneling) of the C2 rotor and the 
uniaxial rotational diffusion of the encapsulating C§4 molecule, the second term oc e~ 3t / T 
accounts for the classical motion of the C2 rotor when its plane of motion is changing 
with the meroaxial reorientations of the encapsulating Cs4 molecule, finally the third term 
oc C(t)e~ 3t ' T accounts for the interference of the two motions of the C§4 molecule with the 
tunneling of the C2 rotor. 

Similarly, using again Eq. (15.81) and the powder-averaged products (A^y) 2 etc. in Ap- 
pendix [HI we find 

(A ZY (t)A ZY (0)) = ^e~*' T , (5.14) 

(B ZY (t)B ZY (0)) = [1 + 2e- 3 ^] , (5.15) 
1 



{C ZY (t)C ZY (0)) = — [ 4 + 17e- 3t/T ] , (5.16) 



and hence 



<«f y (t)af y (0)> = -D(t), (5.17) 

with D(t) again given by Eq. (I5.13p . In the limit of small relaxation time we recover Eq. 
( 14. 12p for static meroaxial disorder. The Raman scattering laws are given by 



Rzzzz(uj) = N(a 2 5(u) + -b 2 D(u)), (5.18) 



and 



b 



2 



Rzyzy(u) = N-D(u). (5.19) 
y 



The Fourier transform of D(t) leads to the scattering function 



3 1 3 

D(u) = C{u) + -J(u) + —G(u). (5.20) 
5 5 



The spectral function C{u) is given by Eqs. (I3.26I) - (I3.3ip while 



J H = ~ ^{l) ^ (5-21) 
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is the Fourier transform of the relaxation function e~ 3t ^ T . The Fourier transform of the 
interference term 

i r+oo 

G(u) = — dte- iu)t C(t)e- m/r (5.22) 
2vr ,/_ 00 

is rewritten as 

/+oo 
dw'C{u-u')J{u'). (5.23) 
-oo 

Using Eqs. (I3.26P - (13.3 1 j) and (15.211) we obtain the scattering function 

G(lo) = G ++ (lu) + G__(u;) + G + _(w) + G_ + (w). (5.24) 

The functions G ++ (c<j), . . . , G (u) have the same structure as C ++ (u), . . . , G (u), Eqs. 

( I3.27P - (I3.30p . respectively, but where the Lorentzians F{uo — w^±i), Eq. (13.311) . are 
replaced by 

r 



(5.25) 



.(^-^ ±1 ) 2 + r2_ 

Similarly to G(o>), Eq. (I3.26p . the function G{uj) is a sum of weighted Lorentzians centered 
around u = ^mrn±\ but °f width 2r where 

T = AD K + 3/r. (5.26) 

The broadening of the transition frequencies of the quantum rotor with increasing temper- 
ature is now due to the uniaxial rotational diffusion and the meroaxial reorientations of 
the encapsulating Cg4 molecule. Notice that both contributions depend on temperature. In 
Fig. [6] we have plotted the function G(u) for several temperatures. The parameters describ- 
ing the dynamics of the Csa molecule are ( n = 100 for the rotational diffusion model and 
Tq 1 = 3 x 10 12 s^ 1 (attempt frequency) and E a = 580 K (activation energy) for the thermally 
activated meroaxial reorientations. Comparable values of the activation energy, i.e. 32(7) 



meV and 35(15) meV have been deduced from neutron scattering studies in solid C70 [17] 



and solid Ceo [271 ] . respectively. While for T = 20 K and 60 K the contribution of 3/r to 
the half width T is negligible in comparison to 4Dr, both (additive) contributions become 
comparable at 150 K. At higher T the thermally-activated reorientations are dominant and 
lead to a smearing out of the low-frequency resonances in the scattering function G{uj). 
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FIG. 7: Spectral function D(u) of the low-frequency Raman scattering laws for T = 20 K, 60 K, 
120 K. 
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If one would artificially exclude the tunneling motion, the function D(uS) entering the 
Raman scattering laws of the C2 unit would reduce to a superposition of Lorentzians centered 
at uj = 0: 

3 1 3 

D(u) = F(u) + -J{lo) + —H{u). (5.27) 
5 5 

The first term on the right-hand side [given by Eq. f)3.24p ] accounts solely for the rotational 
uniaxial diffusion, the second term for the meroaxial reorientations and the last term for the 
interference of these classical motions of the encapsulating C 8 4 molecule. 



VI. DISCUSSION 



It has been shown that the low-frequency (rotational) part of the Raman scattering 
spectrum of a powder crystal of Sc 2 C 2 @C 84 fullerite reflects the superposition of the quantum 
tunneling motion of the encapsulated SC2C2 complex about its long axis and the random 
classical rotational motion of the surrounding C 8 4 molecule. The effect of the C 84 molecule 
on the dynamics of Sc 2 C 2 is twofold. Firstly, since the long axis of Sc 2 C 2 gyroscope coincides 
with the S4 axis of the molecule, the rotation of Sc 2 C 2 about this axis corresponds to the 



motion of the C 2 bond as a planar quantum rotor in a fourfold potential 6]. Secondly, any 
rotation of the C 84 molecule caused by torques due to the thermal lattice environment leads 
to a dragging of the enclosed Sc 2 C 2 unit and hence affects the spectrum of the C 2 quantum 
rotor seen in the laboratory frame. The low-frequency Raman spectra resulting from the 
interaction of the scattering radiation with the induced dipole of the C 2 rotor reflect these 
features. 

In analogy with the dynamics of the C70 molecule in solid C70 , we have assumed that 
the rotational motion of the C 8 4 molecule at a lattice site in a given crystallite is composed 
of two parts: uniaxial rotational diffusion about the S4 axis and stochastic jumps of the 5*4 
orientation among (100) directions. The superposition of the tunneling motion of the planar 
quantum rotor with the classical rotations of the C 8 4 molecule leads to the spectral function 

3 1 3 

D(u) = C(u) + ~J(u) + 4 G (u)> (6- 1 ) 
5 5 

given by Eq. (15.2011 . in the Raman scattering laws Rzzzzi^) and Rzyzy{u). 

The function C(u), defined by Eqs. (13.251) - (13.311) . accounts for tunneling transitions 
between the energy levels of the encapsulated C 2 rotor. The spectrum consists of a series of 
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resonances described by Lorentzians centered at the transition frequencies (Table [Q q = 2) 
and broadened by the uniaxial rotational diffusion (half width 4£)r) of the surrounding Cs4 
molecule. Since the hindering potential for the rotational diffusion about the S4 axis is 
weak, this motion affects the spectrum already at low T. The temperature dependence of 
the spectrum has been studied in Fig. 

The term J(oS) in Eq. (16. ip accounts for the Raman spectrum of the radiation-induced 
C2 dipole while the SC2C2 unit is dragged along by the classical reorientations of the Cs4 
molecule among its three meroaxial directions. This motion which reflects the changes of 
the orientation of the C2 rotor plane is described by a three sites stochastic jump model, 
characterized by a thermally activated relaxation time r = Toe Ea ^ T . Notice that J{oS) leads 
to a central resonance of half width (3/r) in the Raman scattering law even in absence of 
any quantum mechanical tunneling of C2. The width of this central resonance (quasi-elastic 
peak) becomes appreciable at T > 100 K. In the scattering law Rzzzz(u), Eq. (15.181) . this 
quasi-elastic peak is present in addition to the elastic Rayleigh peak. We suggest that in 
future low-energy Raman experiments additional attention will be given to the possible 
identification of the temperature-dependent quasi-elastic peak. 

The last term G(u) in Eq. (16.11) is due to the interference between the uniaxial diffusion- 
modulated tunneling motion described by C{uS) and the stochastic jump model accounted 
for by J(oj). The function G(u) is a convolution of C and J [see Eq. (15.231) ]. While at low 
T the spectra of C(u) and G(u) are very similar (compare the plots for T = 20 K, T = 60 
K in Fig. [5] and Fig. [6]) they become different at higher T (see the 120 K plots) where the 
increasing influence of the stochastic jumps adds to the line broadening. The width 2r of 
the individual resonances, Eq. (15.261) . increases from 1.11 cm -1 at T = 20 K to 3.37 cm -1 
at T = 60 K and 11.43 cm -1 at T = 120 K. This broadening leads to an overlap of the 
low-frequency resonances with increasing T. 

Finally the sum D(u) of these contributions which corresponds to the low- frequency 
Raman response function is shown in Fig. [7J The quasi-elastic peak centered at uj = 
becomes important with increasing temperature. In addition the growing importance of 
G(u>) smears out the low-frequency resonances with increasing T while the higher frequency 
resonances remain prominent. 

The overall shape of the spectral function D(u) and its temperature evolution agree 
very well with the low-frequency Raman scattering results of Ref. [g. There is quantitative 
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agreement with the position of the resonance lines. The smearing out of the low-frequency 
resonances and the prominence of the higher-frequency resonances with increasing T (Fig. 



3 of Ref. 



6) 



are well reproduced by the present theory. In addition to the positions of the 
resonance lines, the theory accounts for their temperature-dependent broadening. In Fig. 
[S]we confront the theoretical spectra D(u) with the experimental Raman spectra, for both 
T = 60 K and T = 120 K. We notice that the experimental spectra are contaminated by a 
plasma line at —uo = 29.6 cm -1 j^]. Note that the central parts of the experimental spectra 
have been omitted in order to remove the effect of the unshifted Rayleigh peak. On the 
other hand the theoretical spectrum exhibits a quasi-elastic peak which is an intrinsic effect 
due to the meroaxial stochastic reorientations of the Sc2C2@Cs4 complex [contribution J{oS) 
in D(u)]. Complementary to the present work it would be useful to measure the dynamics 
of the Cg4 molecule in solid SC2C2C84 directly say by NMR, neutron or /i-spin spectroscopy. 
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APPENDIX A 



We recall that v is the angle of rotation of the Cs4 molecule about its S4 axis. Treating 
the molecule as a Brownian rotor we have the one-dimensional Langevin equation 

QD{t)+cm = m-^- (ai) 

Here B is the moment of inertia about the 5*4 axis, (i> is the friction torque, X(t) a white 
noise driving torque and — dU jdv is a torque due to the orientational crystal field potential 
U{y{t)}. Under the assumption 26] that the friction torque is dominant in comparison to 
;he inertial term and that the variation of U with v is weak, one can use standard methods 



Uj to derive a Smoluchowski equation for the angular distribution function w(y{t)): 



dw ^ d 2 w , . „. 

~dt = (A2) 
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FIG. 8: Comparison of theoretical scattering law D(uS) (solid line), calculated at T = 60 K and 
T = 120 K, with experimental Raman results (dotted line) taken at the respective temperatures. 

where -Dr = k^T/C,. 
The functions 



{^m(^)} 



(A3) 



m = 0, ±1, ... are orthonormal eigenfunctions of the operator d 2 jdv 2 with eigenvalues 
{0, — 1, — m 2 , . . .}. We calculate the conditional probability distribution w{y, t\vo, 0) to find 
the molecule at an angle v at time t when it was at angle Vq at t — 0. The initial condition 
can be written as 



Km w(u,t\u ,0) = 6(u - = y~]ipm( v )A 



(A4) 



where the second member equality is just the closure relation. On the other hand in the 
long-time limit the orientation of the molecule should be random which corresponds to the 
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condition 

lim w(u, t\fo, 0) = — . (A5) 

t >oo 2,71 

A particular solution of (1A2j) subject to these boundary condition is of the form 

w(v,t\v o ,0) = J2^ DRm2t r m (^m^o). (A6) 

m 

The correlation function (cos 2v(t) cos2z/(0)), Eq. fl2T20l) . is rewritten clS cl thermal aver- 
age: 

/*27T /»27T 

F cc (t)= dv du cos2z/G , (z/,t|z/ ,0)cos2z/ . (A7) 
Jo Jo 

Here G(Z/, t|z/ , fydvdvQ is the joint probability of finding the C§4 molecule with orientation 
angle Uq = z/(0) in the interval dv§ initially and orientation v = v{t) in dv at time t. One 
has 

G(u, t\vo, 0) = w(u, t\uo, 0)p(u o ), (A8) 
where the conditional probability w(i',t\i>o,0) is given by 

w (u, t\u , 0) = — V e " m2D R V m( ^ o) , (A9) 

27T * — ' 

m 

m = 0, ±1, . . ., and where p(^o) = 1/27T is the initial equilibrium probability. Carrying out 
the integrals in Eq. (1A7I) gives as a result Eq. (I3.20p . In a similar way we obtain Eq. (13.211) . 

APPENDIX B 

Here we give details about the calculation of the powder averages in Sect. IVl We start 
from the situation where the Csa molecule is in standard orientation, which corresponds 
to the polarizabilities a^ z and a^y given by Eqs. (12. 3p and (12.91) respectively. In order to 
apply the rotation operation R(a,/3, 7) we rewrite the polarizabilities in terms of spherical 



harmonics l^ m (6 l ,0). We use the notations and conventions of Bradley and Cracknell 28]. 
With 

5 n 1/2 



W)=(l6^J (3cos^-l), (Bl) 
Y^ zl (0, <$>) = [ — ) cos6sva.ee***, (B2) 

07T / 
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we get 

1/2 



(!) _ 

a zz ~ a 



+ &(^f) (B3) 
*(^) 1/2[F2l( ^ 0)_F2 " ( ^ 0)]l ^ /2 - (B4) 



1/2 

(1) ' - ' 



The condition <p = tt/2 ensures that the C2 rotor initially moves in the (Y, Z) plane. Al- 
though the function Y^O) does not depend on 0, the condition = 7r/2 has to be taken into 
account after the application of a rotation operation. The transformation law of spherical 
harmonics under a rotation R(a,(3,j) reads: 

+1 

R(a, 0, j)Yr(e, <P)=J2 W 0I>L(a. 0. 7)- (B5) 

n=— I 

Here D^ TO (a,/3,7) are the Wigner rotator functions defined by 

D l nm (a, (3, 7) = C nm e-^d l nm ((3)e- ma , (B6) 



where C nfn = jl*l+nj-M-«» 



281 ] . The functions d l nm (/3) are polynomials in sin(/3/2) and 



cos(/3/2). They satisfy the relations 

dL(0) = (-l) n+m ^-„,- m (/9) = (-l)" +m 4n(/?)- (B7) 

We recall that the angles 9 and on the right-hand side of Eq. (1B5j) refer to the values 
before the application of the rotation. Applying R(a,{3,j) to a^ z , Eq. flB3j) . we have to 
evaluate 

+2 

fl(a,/3, 7 ) K°(0)|^ =7r/2 = E ir(e,7r/2)D^(a,/3 l7 ). (B8) 

n=-2 

In addition to the definitions (1B1I) and (1B2[) we quote 



We further use 



Y± 2 (H.x/2) = - X I^Bin>9. (B9) 



AM = d 2 _ 20 (P) = (76/4) sin 2 /?, (BIO) 
<4(£0 = d -io(P) = -(a/6/4) sin2/3, (Bll) 
4,^) = (3cos 2 /3-l)/2. (B12) 
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We then find after some bookkeeping 

R(a,P, 1 )Y 2 °(9)l =n/2 =(-^-) ' [4 1 l(/?,7) + 4 1 i(A7)sin2^ + C«(/5, 7 )cos2^. 



where 



64tt 



(B13) 



4z(A7) = sin2/5sin 7 , 



cos 2 (3 — sin 2 (3 cos 2 7 — - 



C { zLW, 7) = 2 [ sin2 Z 3 cos 2 7 + 3 cos 2 /3 - 1] . 



(B14) 

(B15) 
(B16) 



Hence 



a 



(1)R 

zz 



+ b[Af z {(3^) + 4z(A7) sin 26 + 7 ) cos 20]. 



(B17) 



With the change (a, /3, 7 ) — >■ (a,/3, 7 — |) we find the coefficients entering 



(2)fl 



4 2 1(A7) = 2 



1 



cos 2 /3 + sin 2 (3 cos 2 7 



4 2 i(A7) = -sin2/5sin 7 , 



4z(A 7) = ^ [" sm 2 /? cos 2 7 + 3 cos 2 /3 - l] 



(B18) 

(B19) 
(B20) 



and with (a, /3, 7 ) — > (a, /3 + |, 7 = 0) the coefficients entering a zz 



(3)B 



4 3 i(^7) = ^ 

>(3) 



sin 2 /3 — cos 2 /3 



£^(/?, 7 ) = 0, 



4 3 l(A7) = ^[cos 2 /3 + 3sin 2 /3-l]. 



(B21) 

(B22) 
(B23) 



The average polarizability a zz is then given by Eq. (14.31) . Applying the definition of powder 

i 



1,2,3; [B% 



average Eq. (14.51) we get [A zz ] 2 

[Czz\ 2 = l|j [Azz Azz} 2 = ~~h-> i 7^ J! [Ezz^i 2 
j =3; [Cf z Cfv 



(0 12 



15 ' 



1,2; [C 



W 12 

zzJ 



15' 1 



1,2, 



_2_ 

15 ' 



In a similar way we apply the rotation operation to a Z y, Eq. (1B4I) . and calculate 



+2 



^0,^)1^(6,^/2) = J2 ^2 n (^,vr/2)^ 1 («,/3, 7 ). 



(B24) 



n=-2 
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Making use of 



d ii(i3) = - cos 2 ^ sin 0, 
d -2i(P) = sin 2 ^ sin/?, 
d 2 m {(3) = v/6/4sin2/3 



d 2 n (f3) = cos 2 | (2 cos /3- 1), 
d 2 n (/3) = sin 2 ^(2cos/3+l), 



(B25) 
(B26) 
(B27) 



and of the relations flB7j) we find 

H)j2(a, /?, 7) [n 1 ^, 0) - nr 1 ^, ^ ,^ /2 

1/2 



(ir\ 1/2 
—J [Ag,(ar, /3, 7 ) + & 7) sin 2£ + C« (a, /3, 7 ) cos 20] (B28) 



where 



. . _ sin 2/5 cos 27 sin a sin 2/3 sin a 
sin p sm 27 cos a H 1 



i?^y (a, /3, 7) = cos /3 cos 7 cos a — cos 2(3 sin 7 sin a, 



<7«(a,/9 )7 ) = 2 



. . „ sin 2/3 cos 27 sin a 3 sin 2/3 sin a 
sm p sm 27 cos a 1 



Hence 



(B29) 
(B30) 
(B31) 



af Y R = b [A ( £(a, /3, 7) + 4r(«> ft 7) + C<g.(a, /3, 7)] • 
With the change of angles (a,/3, 7 — |) we get the coefficients entering afy R : 

4 2 ^(a,/3, 7 ) = 



sin 2/3 cos 27 sin a sin 2/3 sin a 
— sm p sm 27 cos a 1 



? (2) 
'ZY 

cfi(«,A7) 



Bgy^a, (3, 7) = cos (3 sin 7 cos a + cos 2/? cos 7 sin a 
1 



sin 2(3 cos 27 sin a 3 sin 2/3 sin a 
— sin p sin 27 cos a H 1 



and with (at, /3, 7) — > (a;, /3 + |, 7 = 0) the coefficients entering aj^y : 



A y ZY (a, /3,7j = — - sin 2/3 sin a, 

? (3) 



B ( Z y(a, (3, 7) = — sin /3 cos a, 
C^y(a, /3, 7) = — - sin 2/3 sin a. 



(B32) 



(B33) 
(B34) 
(B35) 



(B36) 
(B37) 
(B38) 
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With the definition of powder average Eq. (14. Ill) we get [A^y] 2 — i — 1, 2, 3; [-B^y] 2 = §, 
„• _ i o. [r( 3 )i2 _ i. \r<{i) 12 _ i 7 - _ 1 o. r^(3)i2 _ j_. >iW L -/ i- rW rCj) _ n 

.• ,-. r (i) r w _ _ 1 7 _ 1 o. r (1) r {2) - 1 

4 7^ j; u ZF°zy — 10' 1 ~ z ' zy zy — io - 
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